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ABSTRACT 

We have investigated the steUar content of Barnard 59 (B59), the most active star- forming core in 
the Pipe Nebula. Using the SpeX spectrograph on the NASA Infrared Telescope Facility, we obtained 
moderate resolution, near-infrared (NIR) spectra for 20 candidate Young Stellar Objects (YSOs) 
in B59 and a representative sample of NIR and mid-IR bright sources distributed throughout the 
Pipe. Measuring luminosity and temperature sensitive features in these spectra, we identified likely 
background giant stars and measured each star's spectral type, extinction, and NIR continuum excess. 

To measure B59's age, we place its candidate YSOs in the Hertzsprung-Russell (HR) diagram and 
compare their location to YSOs in several well studied star forming regions, as well as predictions 
of pre-main sequence evolutionary models. We find that B59 is composed of late type (K4-M6) 
low-mass (0.9-0.1 Mq) YSOs whose median stellar age is comparable to, if not slightly older than, 
that of YSOs within the p Oph, Taurus, and Chameleon star forming regions. Deriving absolute 
age estimates from pre-main sequence models computed by D'Antona et al., and accounting only 
for statistical uncertainties, we measure B59's median stellar age to be 2.6±0.8 Myrs. Including 



potential systematic effects increases the error budget for B59's median (DM98) stellar age to 2.6 



+4.1 
-2.6 

Myrs. We also find that the relative age orderings implied by pre-main sequence evolutionary tracks 
depend on the range of stellar masses sampled, as model isochrones possess significantly different mass 
dependences. 

The maximum likelihood median stellar age we measure for B59, and the region's observed gas 
properties, suggest that the B59 dense core has been stable against global collapse for roughly 6 
dynamical timescales, and is actively forming stars with a star formation efficiency per dynamical time 
of ~6%. While the ~150% uncertainties associated with our age measurement propagate directly into 
these derived star formation timescales, the maximum likelihood values nonetheless agree well with 
recent star formation simulations that incorporate various forms of support against collapse, such as 
sub-critical magnetic fields, outflows, and radiative feedback from protostellar heating. 
Subject headings: Hertzsprung-Russell (HR) diagram - infrared: stars - stars: emission-line - stars: 
formation - stars: pre-main sequence 



1. INTRODUCTION 

The current consensus is that star formation pro- 
gresses rapidly over a wide range of environments and 
spatial scales. This view is supported by observa- 
tions that a stellar population's age spread is rarely 
larger th an twice the syst em's crossing time {tctoss = 
Rcore/<yv lElmegreeiil l2000f ) and that ongoing star for- 
mation appears ubiquitous within well known molecu- 
lar clouds: there are few examples of clouds in a pre- 
star forming state, or of clouds associated wit h stel- 
lar populations older than ^5 Myrs (Hartmann et aLl 
l200lUBallesteros-Paredes fc Hartmannll2007| ) . Taken to- 
gether, these observations suggest that star formation 
starts soon after a molecular cloud forms, and persists 
only a few dynamical times before the cloud is disrupted 
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or dissipates. 

The Pipe Neb ula (lAlves e t al. 2007^ and California 
Molecular Cloud (jLada et al.;.2009) are anomalies within 
this prompt star formation paradigm. Lacking ma- 
jor star formation activity, these molecular complexes 
have only recently begun to receive significant scrutiny. 
The first glo bal stu d y of the Pipe was the CO sur- 
vey by lOnishi et al.l ( 1999f ). which identified 14 dense 
clumps that may be fu ture sites of active star formation. 
iLombardi et al.( (I2006D t hen applied the NICER tech- 
nique iLombardi k Alvei (|2001l) to 2MASS data to con- 
struct an extinction map of the Pipe. This map, which 
is shown in Figure 1 , leverages the dense screen of back- 
ground stars provided by the Pipe's projection against 
the Galactic bulge, as well as its proximity (d ^130 pc; 
ILombardi et al.ll2006f ). to resolve physical structures as 
small as 0.04pc. 

Studies of the Pipe's dense cores have begun to in- 
form our understanding of the earliest stages of the 
star formation process. Using a wavelet filtering tech- 
nique to identify individ ual dense cores in t he Lombardi 
et al. extinction map, lAlves et al.l ()2007[ ) noted that 
the Pipe's dense core mass function qualitatively resem- 
bles the stellar mass functions measured in young clus- 
ters ( e.g. iMuench et al 2002; Luhman 200 4) and the 
field (|Covev et all 120081 : [Bastian et al.ll2010[ ). but with 
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Fig. 1. — The locations of MIPS bright sources observed in this program (red diamonds), labehed with the [CLR2010] 
source number f rom T able [2J Provided for context is the NICER extinction map of the Pipe Nebula produced by 
iLombardi at al.l ()2006[ ) (grayscale contours); tha blue box in the upper right identifies tha location of B59, shown in 
mora detail in Figure [H 



a peak mass that is larger by a factor of three. Alvas 
at al. interpreted this as evidence that the stellar mass 
function largely reflects tha physics of core formation, 
with a global star formation efficiency of ^ 33% be- 
tween the dense core and stellar stages. Subsequent 
studies have refined the mass function of Pipe cores, 
adopti ng physically motivate d criteria to identify dense 
cores |Rathborne et al.ll2009[ ). and using detailed tests 
to verify the fidelity and complet eness of the core mass 
function (iKainulainen et all 120091 ). Other authors, how- 
ever, have emphasized the difficulty of identifying co- 
herent physical struc tures from two-dim ensional column- 
density maps (e.g., iSmith et al.l I2008D and cautioned 
against making strict one-to-one co mparisons between 
dense core and stellar m ass functions (jSwift &: Williamj 
[200l[Smith et al.l[2009h . 

The physical properties of Pipe dense cores have also 
shed light on the formation and evoluti on of molecu- 
lar c l ouds. A coordinated radio survev (iMuench et al.l 
120071: iRathborne et al.l 120081: iLada et al.l I2008D identT 
fied thermal pressure as the dominant source of sup- 
port against collapse for the Pipe's dense cores, with 
non-thermal support limited to sub-sonic flows. With 
typical sizes of 0.1 parsec, and measured velocity dis- 
persions of 0.1-0.2 km s~^, these dense cores possess 
sound crossing times and minimum ages of ~1 mil- 
lion years. These core properties, and particularly 
their derived lifetimes, could help discriminate between 



theor e tical models where turbu l ence (iKlessen et all 
[2000t iMac Low fc K lessenI [20 0l lOffner et all 120081) 
or magnetic fields (Naka nol 119841: IShu et'aLl Il987t 
iMouschovias fc Ciolekiil999t iNakamura fc Lil 120081 ) pro- 
vide the dominant means of support for molecular cores. 
This evidence, however, is somewhat circumstantial, and 
relies on a critical assumption: that the Pipe's dense core 
population will one day become active sites of ongoing 
star formation. 

The Barnard 59 dense core (hereafter B59 : originally 
dubbed the 'Sinkhole' by iBarnard et al.l[l927l ) could pro- 
vide a clearer test for these theoretical models: as the 
only site of active star formation within the Pipe, no 
extrapolations of the core's future behavior are neces- 
sary. While B59 is the Pipe's most massive core, its 
physical properties (including the dominance of thermal 
motions) are nonetheless representative of the broader 
population of dense cores within the Pipe, making it an 
instructive case study. Since 1950, seven optically visi- 
ble Ha emitting stars (or unresolved multiple systems) 
have been associated with B59 (V359 Oph, KK Oph AB, 
LkHa 345 Oph, LkH a 346 Oph AB, KW 002 THa 27-5 
and LkHa 347 Oph : Merrfll & Burwell 1950; The 196^ 
[Stephenson fc Sanduleak 1977 : Kohoutek &: Wchmeve: 
120031 iHerbid I2005D. K recent Spitzer survey of B59 
bv iBrooke et al.l (|2007[ ) identified 16 additional candi- 
date members with strong near-infrared (NIR) excesses, 
4 of which were detected previously in the sub-mm 
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(jReipurth et al.|[T996D or by IRAS at far infrared wave- 
lengths. These sources appear to span a range of evo- 
lutionary states, from the low -mass, relat i vely e mbedded 
Class O/I protostar studied bv lRiaz et al.l ()2009( ). to more 
revealed Class II sources. Analyses of Spitzer, ROSAT, 
and XMM observations of the Pipe Nebula confirm that 
B59 contains the vast majority of the Pipe's recently 
formed stars (|Forbrich et aL|[200l [20Tol ). 

The physical properties of the B59 dense core make it a 
good analog for modern computational simulations of the 
collapse of isolated dense cores. In this respect, charac- 
terizing the timescale and efficiency of star formation in 
the B59 dense core would provide a good opportunity to 
compare these empirically measured quantities with the 
predictions of these simulations. Various methods have 
been developed to estimate the ages of pre- main sequence 
stars, with differing advantages and levels of accuracy. 
The destruction of a star's primordial lithium abundance 
has been used to estimate the ages of pre-main sequence 
stars with moderate success (i.e., typical uncertainties 
of ~10 Myrs for me mbers of 10-30 Myr moving groups; 
iMentuch et al.ll2008[ ). but this method would be difficuh 
to apply to B59's YSOs: measuring lithium abundances 
requires high resolution, high S/N optical spectra that 
would be difhcu h to acquire for B5 9's heavily extincted 
YSOs (Ay > 5 lBrookeeraIll2007[ ). and recent studies 
indicate that lithium depletion may not a good a ge indi- 
cator for young (t<10 Myrs) T Tauri stars (Sestit o et al.l 
|2008^. A pre-main sequence star's surface gravity also 
rises as it contracts towards the main sequence, provid- 
ing an additional opportunity to constrain the star's age. 
Theoretical models suggest that a precision of 0.1 dex in 
log g can provi de age estimates accura te to within 50% 
(see Figure 8, iDoppmann et al.l 120051 ) . but this preci- 
sion is difficult to achieve in practice. Log g values have 
been measured for pre-main sequence stars with an ac- 
curac y of ~0.2 dex ([Doppmann et al.ll2005t iTakagi et al.l 
120101 e.g., ), providing age estimates accurate to ~l-2 
Myrs, but as with lithium depletion studies, these mea- 
surements depend on precise measurements of line shapes 
and strengths, requiring observationally intensive high 
resolution, high signal-to-noise spectra: obtaining spec- 
tra of this quality for ^20 B59 YSOs would require mul- 
tiple nights on an 8-meter class telescope. Surface grav- 
ity indicators have been developed for lower-resolution 
spectra, but these indicators typically depend on obser- 
vations at J or shorter wavelengths, again difficult to ob- 
tain for B59 YSOs, and have a typical precision of ±0.5 
dex, providing somewhat poorer age discrimination than 
is ava .ilable with standard HR diagram ana lyses {^ 3-5 
Mvrs lGoriova et al.l [20031: 1 AUers et all 120091) . 

To characterize the timescale, efficiency, and products 
of star formation within the Pipe Nebula, we have con- 
ducted the first spectroscopic census of B59's stellar pop- 
ulation. Our observations discriminate bona-fide B59 
members from background giants in the Galactic bulge, 
and establish that these members are late type (K4-M6) 
low-mass (0.9-0.1 Mq) young stellar objects (YSOs). We 
then derive estimates of the absolute and relative ages of 
B59's YSOs by comparing their locations within the ob- 
servational Hertzsprung-Russell diagram with theoreti- 
cal models and YSOs in other well known clusters. We 
measure B59's median stellar age to be ~2.6l2'6 Myrs; 
while the error bars on this measurement are significant. 
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Fig. 2. — The K band image of B59 obtained by 
iRoman-Zuniga et al.l ()2009l ). with locations of candi- 
date YSOs observed in this program iden tified by dia- 
monds and the source number assigned bv lBrooke et al.l 
(2 007). The NICER-ba s cd ex tinction map constructed 
by IRoman-Zuhiga et al.l (|2009f ) is shown for comparison 
as white contours. The lowest contour represents an ex- 
tinction of Ay =2.5; contours correspond to a 5 mag. 
increase in extinction up to Ay =22.5, and increases of 
10 mag. thereafter. 

the maximum likelihood estimate implies that star for- 
mation has persisted within B59 for more than five times 
the current dynamical timescale and crossing time of 
B59's dense gas. 

In Section 2, we present our SpeX observations of can- 
didate B59 members. Our method for measuring stellar 
parameters from these spectra is introduced in Section 3, 
along with our algorithm for placing those objects in the 
HR diagram and inferring ages and masses from pre-main 
sequence evolutionary models. We discuss our results in 
Section 4, and summarize our findings in Section 5. We 
also include a NIR spectral atlas of TW Hydrae members 
as an appendix to this work. 

2. OBSERVATIONS 

Using the SpeX spectrograph (jRavner et al.ll200"3l ) at 
NASA's 2.5 meter Infrared Telescope Facility, we ob- 
tained low (R~200) and moderate (R ~ 2000) resolu- 
tion spectra of candidate young stellar objects (YSOs) 
in B59 and the Pipe Nebula. Observations were carried 
out during the early mornings of April 14-18 2008 (HST); 
typical seeing was 0.7", with clear skies interrupted by 
periods of light cirrus. The bulk of our observations were 
carried out in SpeX's SXD mode with the 0.3"slit, pro- 
viding nearly contiguous spectral coverage from 0.8 to 2.5 
microns with a resolution of R~2000. During a night of 
abnormally poor seeing (~1") we made use of the larger 
0.5"slit, degrading the spectral resolution to R~1200. 



Covey et al. 



TABLE 1 
IRTF/SpeX SXD Observing Log for Candidate B59 Members 



Source 






2MASS 


2MASS 


2MASS 


SXD 


SXD 


SXD 


SXD 


SXD 


NTT 


NTT 


ID 


RA 


Dec 


J 


H 


Ks 


itime 


slit 


J S/N 


H S/N 


KS/N 


H 


K 




BHB2007 


1 


17:11:03.93 


-27:22:55.08 


10.46 


8.99 


7.76 


240 


0.3 


152/52 


243/49 


293/114 


9.29 


8.29 




BHB2007 


2 


17:11:04.12 


-27:22:59.44 


9.75 


8.79 


8.05 


240 


0.3 


136/48 


155/35 


151/73 


8.91 


8.56 




BHB2007 


3 


17:11:11.83 


-27:26:54.98 


14.09 


12.62 


11.76 


6240 


0.5/0.3 


58/13 


108/26 


186/40 


12.69 


11.77 




BHB2007 


4 


17:11:14.46 


-27:26:54.45 


11.62 


10.42 


9.63 


1920 


0.3 


277/41 


344/48 


347/57 


10.62 


9.77 




BHB2007 


5 


17:11:16.26 


-27:20:28.82 


11.23 


9.70 


9.07 


1440 


0.3 


132/57 


168/23 


158/56 


9.70 


9.08 




BHB2007 


6 


17:11:16.31 


-27:25:14.50 


10.58 


9.37 


8.75 


960 


0.3 


203/40 


234/41 


227/63 


9.37 


8.75 




BHB2007 


7 


17:11:17.28 


-27:25:08.16 


13.61 


10.82 


8.77 


1920 


0.3 


58/12 


261/34 


325/92 


10.58 


8.77 




BHB2007 


8 NW 


17:11:18.24 


-27:25:48.94 








1920 


0.3 




2/0 


59/17 


16.01 


12.87 




BHB2007 


8 SE 


17:11:18.32 


-27:25:49.53 


>18.4 


15.12 


11.95 


1920 


0.3 




3/0 


60/16 


15.35 


12.78 




BHB2007 


9 


17:11:21.55 


-27:27:41.70 


12.74 


10.57 


8.98 


2400 


0.3 


122/31 


230/42 


348/69 


11.215 


8.98 




BHB2007 


12 


17:11:25.08 


-27:24:42.67 


16.52 


13.61 


11.57 


2880 


0.5 


15/3 


63/13 


154/34 


13.32 


11.58 




BHB2007 


13 NW 


17:11:26.95 


-27:23:48.40 


11.88 


10.14 


9.08 


1440 


0.3 


84/32 


126/42 


120/3 


10.14 


9.08 




BHB2007 


14 


17:11:27.28 


-27:25:28.32 


13.17 


10.65 


9.11 


2400 


0.5 


120/21 


351/31 


494/111 


10.57 


11.44 




BHB2007 


15 


17:11:29.43 


-27:25:36.89 


13.28 


11.74 


10.70 


2400 


0.5 


149/24 


291/30 


359/46 


11.74 


10.58 




BHB2007 


16 


17:11:30.35 


-27:26:28.91 


11.91 


10.00 


8.89 


1680 


0.3 


49/19 


64/32 


64/66 


10.23 


9.19 




BHB2007 


17 


17:11:41.00 


-27:18:36.82 


10.64 


8.75 


7.54 


960 


0.3 


146/18 


150/8 


106/19 


8.75 


7.54 




BHB2007 


18 NE 


17:11:41.73 


-27:25:50.26 








480 


0.3 


216/37 


287/39 


312/50 


10.46 


9.86 




BHB2007 


18 SW 












1440 


0.3 


72/32 


76/40 


69/55 








BHB2007 


19 


17:11:43.17 


-27:30:58.54 


14.09 


12.75 


11.91 


1920 


0.3 


47/11 


92/17 


127/26 


12.66 


11.74 




BHB2007 


20 


17:12:00.20 


-27:20:18.10 


10.54 


9.64 


9.17 


960 


0.3 


119/51 


116/51 


146/81 


9.641 


9.17 


V359 Oph 


17:08:54.28 


-27:12:33.03 


9.73 


9.12 


8.80 
















LkHa 345 


17:10:48.04 


-27:40:51.19 


9.60 


8.69 


8.35 
















KK Oph A 


17:10:08.07 


-27:15:18.24 


9.07 


7.23 


5.80 

















For maximum observational effieiency, we chose not 
use the image rotator to aUgn the sht with the parallac- 
tic angle over the course of each night. Spectra of telluric 
standards taken at slit rotations that differ by 30-60° sug- 
gest that color-dependent slit losses due to atmospheric 
refraction are relatively minor, altering the slope of the 
spectral energy distribution (SED) between the J and K 
bands by less than 10%. During our last night of ob- 
servations we also utilized SpeX's high-throughput, low- 
resolution (R^200) prism mode to obtain high fidelity 
NIR SEDs for candidate YSOs in B59. During these ob- 
servations, care was taken to align the slit to within 10° 
of the parallactic angle, except in cases where the slit was 
rotated to allow observations of individual components of 
a n apparent multiple system. 

iBrooke et aP (|2007[ ) identified candidate YSOs in B59 
by searching for IR excesses in spectral energy distribu- 
tions (SEDs) constructed from 2MASS, IRAC and MIPS 
photometry. For the remainder of this paper, we iden- 
tify these candidate YSOs with a [BHB2007] prefix, fol- 
lowed by the source ID assigned by Brooke et al. in their 
Table 1. We obtained spectra for each of these candi- 
date B59 members, with the exception of two objects 
too faint to observe with SpeX ([B HB2007] 10, fc 11). A 
K ban d image of B59 observed by iR oman-Zuniga et al.l 
()2009( ) is shown in Figure [2l overplotted with the loca- 
tions of the candidate YSOs observed in this program. 
Our SXD observations are summarized in Table [U where 
we indicate the slit width and exposure time used to 
obtain each source's spectrum. We include for com- 
pleteness the 2MASS position and magnitudes for each 
source as well as H and K b and magnitudes measured by 
iRoman-Zuniga et al.l (|2OO90 from their NTT images of 
the cluster. We also include positions and 2MASS mag- 
nitudes for 3 likely B59 members (V359 Oph, KK Oph 
AB, and LkHa 345 Oph) with spectral types in the liter- 
ature that suffice to place each star in the HR diagram. 

To characterize the more extended population of bright 



IR objects that extends across the Pipe, we obtained 
spectra for a sample of objects with K<10 (to ensure 
accessibility to SpeX) and detections at both 24 and 70 
microns in an extensi ve Spitzer survey of the Pipe Nebula 
(jForbrich et al.ll2009| ) . The locations of these sources are 
indicated in Figure [I] and their NIR properties are listed 
in the top tier of Table [5] We observed all four sources 
([CLR2O1O0 1, 2, 5 and 7) that met our flux criteria and 
are located within high extinction peaks identified by 
lAlves etall (|2007| ). as we considered these sources the 
most likely bona fide YSO candidates. We additionally 
observed 6 similar sources that are not associated with 
extinction cores ([CLR2010 ] 3, 4, 6 8 and 9). Accord- 
ing to the criteria of iForbrich et al.l (|2009f ). two of the 
sources in our sample ([CLR2010] 1 & 2) are MIPS ex- 
cess sources. A third source, [CLR2010] 7, has a MIPS 
spectral i ndex marginally qualif ying as class II but is not 
listed by IForbrich et al.l ()2009D due to a position mis- 
match at the source detectio n stage. 

The SpeX spectral atlasQ fiCushing et al.ll2005l Rayner 
et al., in prep; Vacca et al., in prep.) samples many 
spectral types and luminosity classes. To supplement this 
grid of standards with stars possessing pre- main sequence 
surface gravities, we observed members of the TW Hy- 
drae moving group and p Ophiuchus cluster. The pre- 
main sequence standards are listed in the bottom tier of 
Table [2l and spectra of the TW Hya members are shown 
in Appendix El 

All data were reduced using SpeX's dedicated I DL re- 
duction package, SpeXtool (G ushing et alJ [20M ) . Tel- 
luric absorption was removed using observations of 
AOV stars at airmasses similar to those of the target 

^ For consistency with the [BHB2007] designation adopted for 
the candidate B59 members identified by Brook e et al] 1)20071 ). we 
adopt a source designation of [CLR2010], followed by the source 
number given in Table (2] for the IR bright sources we observed 
throughout the Pipe. 

^ available online at 

http://irtfweb. if a.hawall .edu/$\slm$spex/spexllbrary/IRTFllbrary .html| 
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TABLE 2 
IRTF/SpeX SXD Observing Log for Pipe/TW Hya Sources 



Source 






2MASS 


2MASS 


2MASS 


SXD 


SXD 


SXD 


SXD 


SXD 


ID 


RA 


Dec 


J 


H 


Ks 


itimc 


slit 


J S/N 


HS/N 


KS/N 




CLR2010 


1 


17:14:51.08 


-27:22:39.96 


10.26 


8.89 


8.16 


360 


0.3 


447/55 


654/54 


498/84 




CLR2010 


2 NE 












480 


0.3 


100/44 


107/33 


106/36 




CLR2010 


2 SW 


17:14:56.49 


-27:31:08.08 


10.42 


9.56 


9.04 


480 


0.3 


86/32 


100/32 


108/35 




CLR2010 


3 


17:15:59.37 


-27:00:24.75 


6.99 


5.72 


5.13 


120 


0.3 


507/116 


794/19 


935/58 




CLR2010 


4 


17:19:06.05 


-27:06:18.95 


13.07 


9.69 


7.50 


360 


0.3 


25/5 


136/9 


141/38 




CLR2010 


5 


17:20:12.47 


-26:55:09.40 


9.12 


7.58 


6.67 


720 


0.3 


229/30 


241/10 


294/24 




CLR2010 


6 


17:22:31.01 


-23:53:48.76 


7.52 


6.06 


5.13 


120 


0.3 


357/17 


474/9 


884/21 




CLR2010 


7 


17:26:15.14 


-26:44:37.31 


9.16 


7.13 


5.79 


240 


0.3 


180/64 


180/20 


182/62 




CLR2010 


8 


17:27:41.56 


-26:53:46.7 


10.12 


7.87 


6.19 


240 


0.3 


38/15 


146/43 


488/49 




CLR2010 


9 


17:29:54.08 


-25:59:12.79 


7.29 


5.83 


5.05 


120 


0.3 


423/78 


726/18 


971/56 


TW Hya ] 


2 


11:21:05.6 


-38:45:16 


9.00 


8.33 


8.05 


480 


0.3 


343/113 


392/46 


346/88 


CD 29-8887 A 


11:09:14.0 


-30:01:39 


7.85 




7.18 
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-41:36:39 
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0.3 


227/99 
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16:26:23.7 
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0.3 
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435/57 
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0.3 
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0.3 
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-37:31:53 
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11:10:28.05 


-37:31:54.2 
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(Asecz ^0.1). Hydrogen features were removed from the 
AOV spectra by dividing by a model of Vega, and the 
target's true spectral slope and absolute flux densities 
were recovered by multiplying the target spectrum by a 
~10,000 K blackbody spectrum scaled t o match that of 
the telluric standard (jVacca et al.l 120031 ) . Tables [1] and 
[2] list the signal-to-noise (S/N) levels measured in the 
J, H and Ks bands of each source's reduced SXD spec- 
trum. Two S/N estimates are listed for each band: the 
first estimate is the S/N calculated directly by SpeXtool. 
The second estimate is an empirical lower limit to each 
spectrum's S/N generated from the standard deviation of 
relatively featureless regions in each band (J: 1.215-1.235 
^m, H: 1.52-1.54 //m, K: 2.13-2.15 ixm). 

3. ANALYSIS 

We have analyzed these SpeX spectra to distinguish 
bona fide B59 members from contaminating background 
giants, and to infer the physical properties of each YSO. 
In §3.11 13.21 and 13.31 we describe an iterative technique 
that identifies the spectral type, extinction and veiling 
that best reproduce the individual features and over- 
all shape of each candidate B59 YSO's spectrum. In 
i j3.41 we supplement these spectral measurements with 
photometry from the literature to estimate the effective 
temperature (Tg/^) and bolometric luminosity {Lboi) of 
each target. Placing these YSOs in an observational 
Hertzsprung-Russell (H-R) diagram, we used pre-main 
sequence evolutionary models to infer masses and ages 
for each object. We present in §3.51 an estimate of B59's 
absolute age, as well as a relative age estimate derived 
from comparison to other clusters with similarly homo- 



geneous photometric and spectroscopic datasets. 

3.1. Estimating Extinction & Veiling From Synthetic 
Photometry 

Following fMever et all ()1997[) . we assume that a YSO's 
observed NIR color (e.g., (J — H)obs) is due solely to its 
intrinsic photospheric color (e.g., {J — H)int), an addi- 
tional continuum excess from its circumstellar disk (e.g., 
A(J - H), such that (J - H)cTTS = [J - H),nt_ + 
A(J — H)), and an independent extinction/reddening 
term (e.g., E{J - H) = 0.613 x Kr] [J - H)obs = 
{J — H)cTTS + E(J — H)). Algebraically combining 
the iMever et all (|1997D CTTS locus (transformed onto 
the 2MASS photometric system) with a well behaved 
extinction la\\0 then requires a unique combination of 
reddening and disk excesses/veiling to reconcile a YSO's 
intrinsic and observed (J-H) and (H-Kg) colors: 



Ah = 2.512 X (J-H)obs~ 
1.527 X {H -Ks)obs- 1.254 

A{J - H) = [J - HUs - 
0.613 X Ah - {J - H),nt 



(1) 



(2) 



^ IHoman-Zuniga et a l. (2007) measured the near/mid-IR extinc- 
tion law in B59 itself, identifying Aj/Ak and Ah /Ak ratios that 
agree within Icr of the values measured bv llnd ebetouw et al. ( 200^) 
from large numbers of 2MASS stars. Given th e agreement between 
these two studies, we adopt the llndebetouw "et al.l (2005) extinction 
law for maximum compatibility with the 2MASS photometry we 
utilize to characterize YSOs in B59 and other star formation re- 
gions. 
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TABLE 3 

Spectral Line Indices 



Line 
Name 



Line 
Center {/i) 



Line 
Width ifi) 



Cont. 1 
Center (fi) 



Cont. 1 
Width ifi) 



Cont. 2 
Center (jj.) 



Cont. 2 
Width (At) 



Mg 1.50 


1.504 


0.004 


1.498 


0.0045 


1.5095 


0.0045 


K L52 


1.5172 


0.004 


1.5105 


0.004 


1.523 


0.004 


Mg 1.58 


1.576 


0.004 


1.57 


0.004 


1.58 


0.004 


Si 1.59 


1.59 


0.005 


1.586 


0.003 


1.594 


0.003 


CO 1.62 


1.61975 


0.0045 


1.615 


0.003 


1.628 


0.003 


CO 1.66 


1.6625 


0.004 


1.6565 


0.003 


1.67 


0.003 


Al 1.67 


1.674 


0.006 


1.659 


0.003 


1.678 


0.002 


1.70 dip 


1.7075 


0.003 


1.704 


0.003 


1.7145 


0.003 


Mg 1.71 


1.7115 


0.003 


1.704 


0.003 


1.7145 


0.003 


Ca 1.98 


1.982 


0.013 


1.9676 


0.005 


1.99775 


0.005 


Na 2.21 


2.2075 


0.007 


2.195 


0.004 


2.217 


0.006 


Ca 2.26 


2.264 


0.007 


2.258 


0.0045 


2.27 


0.0045 


CO 2.3 


2.30375 


0.0225 


2.288 


0.007 


2.3185 


0.004 


CO 2.34 


2.3455 


0.003 


2.3425 


0.003 


2.349 


0.003 

















0.354 xAh-{H- Ksh^t 



(3) 



YSOs often demonstrate considerable photometric 
varia bility, including in the NIR regime we focus on here 
(e.g., lCarpente r2001: Moralcs-Caldcron ct al. 2009). To 
estimate each YSO's extinction and NIR excess at the ex- 
act epoch that the spectra were obtained, we calculate 
synthetic NIR photometry from each target' s SpeX spec- 
trum f ollowing the prescription presented bv lCovev et al.l 
(|2007f ) in their §3.3.1. Using these JHK magnitudes, and 
assuming intrinsic photospheric colors characteristic of 
K7/M0 type stars, we can therefore obtain a good initial 
estimate of each YSO's extinction (A//) and A(J-II) and 
A(H-Ks) color excesses. 

3.2. Extinction and Veiling Resistant Spectral 
Classification 

Numerous investigators have identified H and K band 
spectral indices sensitiv e to a star ' s lum i nosity class and 
spectral ty pe (Kleinmann & Hall 1986: Ali ct al. 1995: 
Wallace fc Hinklc 1997: Mevcr ct al. 1998: Aspin 2003: 
Ivanov et al.ll2004HAllers et al.ll2007f )'. We present in Ta- 
ble [3] a set of indices well suited for measuring stellar 
spectral types from moderate resolution (R^lOOO-2000) 
spectra. Each feature's local continuum is measured from 
a linear fit to the median flux in two neaby regions ('Cont. 
1' and 'Cont. 2' in Table E]). Integrating the differ- 
ence between the expected continuum and the actual 
flux in the line (defined by the 'Line Center and 'Line 
Width' columns of Table |3l), and dividing by the contin- 
uum flux density in the center of the line region, provides 
an Equivalent Width (EqW) measure of the strength of 
each feature. Equivalent width measures are convenient 
for spectral classification of YSOs: equivalent widths are 
relatively insensitive to extinction, which affects both the 
line and the local continuum equally, and ratios of closely 
separated lines are relatively robust to excess continuum 
flux, as both lines will be diluted equally and the sum 
effect will cancel out. 

The strength and shape of water absorption in the 
H and K bands has also been identified by numerous 
authors as a useful spectral type and luminosity class 
indicator (e.g., IWilking et al.. ,1999, : ,Rcid et al.. ,2001; 



Geballe et al.l 120021: iMcLean et al.l 120031: iSlesnick et al.l 
2004 lAllers et al.l l2007t [Weights et al.ll2009[) . but is less 
amenable to characterization as an Equivalent Width. 
We defined two indices to characterize water absorption 
in the H and K bands; these indices are calculated as: 



Index = 



< Bandl > / < Band2 > 

< Band2 > / < Band3 > 



(4) 



where <X> simply denotes the mean flux level within 
the range of each band as given in Table |31 

To account for any residual sensitivity to extinction or 
veiling in our spectral indices, we compared our target's 
index values to those measured from artificially reddened 
and veiled spectral standards. As described in ^ 33.11 syn- 
thetic photometry from each SpeX spectrum produces 
an initial estimate of each target's extinction and (J- 
II)/(II-K) color excesses assuming the standard K7/M0 
spectral type. Adopting the characteristic J band veiling 
(i'j=-^ = 0.25) measured bv iCieza et al.l (|2005) for a 
large sample of CTTSs sets the initial absolute scale for 
the veiling fluxes implied by each source's A(J-II) and 
A(II-Ks) color excesses: 



A J = 2.5 X Log{rj + 1) 
Ai/ = AJ + A{J-H) 

AH 

rn = 10^^ - 1 



(5) 



We then iteratively adjusted each target's assumed ab- 
solute J band veiling and spectral type (and thus intrinsic 
photospheric colors) to identify the self-consistent spec- 
tral type/extinction/veiling combination that best repli- 
cated the target's measured spectral indices and overall 
SED shape. Figures [3] and [4] demonstrate our two dimen- 
sional spectral classification of GY 292 and [BHB2007] 
5, respectively. GY 292 appears to be moderately veiled 
(A(J-H): 0.22; A(H-K,): 0.35), but [BHB2007] 5 shows 
no evidence for veiling, so these two stars demonstrate 
the sensitivity of these spectral indices to the presence of 
veiling. Veiling has the largest impact on the two water 
indices shown in the upper left panels of Figures [3] and 
01 standards veiled to match GY 292 are significantly 
compressed compared to the unveiled standards used to 
type [BHB2007] 5. 
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Fig. 3. — Two dimensional spectral classification of GY 292: each panel shows luminosity and gravity sensitive H 
and K band spectral indices presented in Wd.2 Ind ices measured from GY 292, a confirmed YSO with a previously 
measured K8 spectral type (|Luhman fc Riekelll999D . are shown as a filled black point, while indices measured from 
spectral standards are labeled according spectral type, and color-coded according to luminosity class. Standard spectra 
were artificially veiled and extincted to match estimates from GY 292 's photometry and spectra before indices were 
measured. GY 292 appears to be an ~K5 object with subgiant-like surface gravity, as expected for a bona-fide YSO. 



TABLE 4 
Water Indices 



Name 



Band 1 



Band 2 



Band 3 



H2O-H 1.595-1.615 
H2O-K 2.18-2.20 



1.68-1.70 1.76-1.78 
2.27-2.29 2.36-2.38 



Figure S] also demonstrates the utility of two- 
dimensional spectral classification for identifying back- 
ground giants, as [BHB2007] 5 is clearly identified in 
each plot as a giant star; lacking significant veiling, or 
emission lines indicative of mass accretion, this object 
appears to be a background late K/early M giant star. 

3.3. Derived Spectral Types 

Applying our classification algorithm to these spectra 
identifies two of the B59 candidate YSOs and the major- 
ity of the IR excess stars throughout the Pipe as likely 
giants: we classify [BHB2007] 5, [CLR2010] 3 and 9 as 
reddened giants or supergiants, [CLR2010] 8 as a giant 
carbon star, and [CLR2010] 4, 5, 6 and [BHB2007] 17 as 
OH/IR stars, likely residing in the Galactic Bulge. The 
spectra of these stars are shown in Figure [5] along with 



spectr a of known giants obtained by iLancon fc Wood! 
dioOd) and iCushing et al.l (J2005). Our classification of 
[BHB2007] 5 and 17 as background giants is reinforced 
lay their non-detection at X-ray wavelengths in an XMM 
pointing containing B59 (Forbrich et al. 2010, in press). 
The spectra of the confirmed B59 YSOs are shown in 
Figures [6] and [3 with their derived spectral types, extinc- 
tions, and JHK veilings presented in Table [S] We tested 
the accuracy of these results by comparing the parame- 
ter s we measure for four p Oph YSOs to those reported 
bv lLuhman fc Rieki ()1999D . We derived spectral types 
for G Y 292 (K5), GY 17 (K5) GY 250 (Ml) and GY 93 
fM5): lLuhman fc Riekd ()1999f ) derived spectral types for 
these objects of KB, K6, MO, and M4 respectively. The 
types assigned to GY 17, 93, and 250 agree within 1 sub- 
class, while the type we assign to GY 292 is earlier by 
two/three subclasses (following iKirkpatrick et al.l I1991L 
our classification scheme omits types K6 and KB). On 
the basis of this comparison, we adopt a conservative un- 
certainty of ±2 subclasses in our derived spectral types. 
Similarly, comparing the e xtinctions and veihngs w e de- 
rive with those reported by lLuhman fc Rieki (|1999f ) sug- 



gest that we achieve (tah 
Two sources, [CLR2010] 



0.5 mag and ar^ ~ 0.35. 
1 and 7, present somewhat 
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Fig. 4. — The two dimensional spectral classification plots for candidate YSO [BHB2007] 5; symbols are as in Figure 
[3l [BHB2007] 5 appears to be a mid-late K star, with giant-like surface gravity; visual inspection of its spectrum 
confirms it as a background giant rather than a bona- fide YSO (see Fig. [5]). 



ambiguous spectra. [CLR2010] 7 possesses a sub- 
giant/giant-like surface gravity, with prominent Paschen 
/3, Paschen 7, and Brackett 7 emission features as seen in 
some YSOs. [CLR2010] 7 is also located on the edge of 
a Pipe extinction c ore, and is on the lower boundary de- 
fined bv iForbrich et al.l ()2009 !) to select candidate class 
II sources based on their MIPS spectral index. While 
these properties a re consistent with a YSO classification, 
lOiha et al.l (pOOl identify [CLR2010] 7 as a likely high 
mass-loss AGB star in the Galactic Bulge in their anal- 
ysis of detections by the Mid-course Space Experiment. 
This classification is also consistent with the observed 
properties of [CLR2010] 7, explaining its strong HI emis- 
sion as a signature of strong stellar winds. We tenta- 
tively classify this source as a giant star, and include its 
spectrum in Figure [3 [CLR2010] 1 is more ambiguous: 
this source possesses a low (sub-giant) surface gravity, 
and has a location in the HR diagram consistent with 
that of other B59 members (see §3.3 for a discussion of 
HR diagram placement). This HR diagram location im- 
plies a rather young age for [CLR2010] 1, however, which 
is somewhat inconsistent with the source's lack of NIR 
excess/disk emission or mass accretion sensitive emission 
lines. Indeed, [CLR2010] I's NIR photometry is also con- 
sistent with that of a reddened background giant/sub- 
giant. Definitively resolving [CLR2010] I's status will 
require additional observations; for now we include it in 



Table [5] and Fig. |6]as a potential Pipe YSO, but exclude 
it from our population analysis of B59 members. 

One of the IR excess stars located outside of B59, how- 
ever, is likely a newly identified YSO. [CLR2qiO] 2 lies 
in the stem of the Pipe, within an IAlves et al.l ()2007| ) ex- 
tinction core (see Fig. [J for i ts exa ct location) , and was 
identified by IForbrich et al.l ()2009[ ) as a new candidate 
YSO (source 16 in their object list). While detected as a 
single source by MIPS, [CLR2010] 2 proved to be a visual 
binary when observed with SpeX: we obtained spectra 
of both the NE and SW components, and the spectral 
indices measured from both sources are consistent with 
pre-main sequence surface gravities. Strong HI, Hel, and 
Ca II emission lines are also visible in the spectrum of 
[CLR2010] 2 SW, consistent with the source being an 
actively accreting YSO. This source was also identified 
as an Ha emission line source by Kohoutck & Wehmeyea 
(2003) , indicating that the HI emission persists over long 
timescales. These observations suggest that [CLR2010] 2 
SW is a bona fide YSO associated with the Pipe Nebula, 
and [CLR2010] 2 NE is a promising candidate as well. 
The detection of Ha from [CLR2010] 2, and the low ex- 
tinction estimates we derive from our SpeX observations, 
indicate that follow-up optical spectroscopy of the 6 70 7 A 
Li line may be a tractable means of definitively resolv- 
ing each component's pre-main sequence status. As with 
[CLR2010] 1, we include [CLR2010] 2 NE and SW in 
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TABLE 5 
Measured Spectral Parameters 



Source 



SpT 



rj 



AJ-H TH AH-K XK Br7? 



low 
log g? 



YSO 



BHB2007 
BHB2007' 
BHB2007' 
BHB2007' 
BHB2007' 
BHB2007 
BHB2007 
BHB2007 
BHB2007 
BHB2007 
BHB2007 
BHB2007 
BHB2007 
BHB2007 
BHB2007 
BHB2007 
BHB2007 
BHB2007 
CLR2010' 
CLR2010' 
CLR2010^ 



1 


K7 


0.84 


0.26 


0.05 


0.32 


0.30 


0.76 


Y 


Y 


Y 


2 


MS 


0.68 


0.00 


0.00 


0.00 


0.00 


0.00 


W 


Y 


Y 


3 


M6 


1.43 


0.00 


0.00 


0.00 


0.00 


0.00 


7 


N 


7 


4 


M5 


0.74 


0.26 


0.10 


0.39 


0.16 


0.62 


Y 


Y 


Y 


6 


M2 


0.97 


0.00 


0.00 


0.00 


0.00 


0.00 


N 


Y 


Y 


7 


K5: 


2.80 


0.45 


0.44 


1.16 


0.69 


3.11 


Y 


Y 


Y 


8NW'' 


K/M: 


2.41 






1.20 


1.54 


11.5 


N 


7 


7 


8SEb 


M: 


4.54 






0.26 


3.89 


1.34 


N 


7 


7 


9 


K5 


1.S4 


0.38 


0.39 


0.97 


0.61 


2.48 


Y 


Y 


Y 


12 


M5 


3.00 


0.26 


0.13 


0.42 


0.27 


0.73 


Y 


Y 


Y 


13NW 


M2 


2.18 


0.15 


0.00 


0.15 


0.00 


0.15 


W? 


Y 


Y 


14 


K5 


2.48 


0.38 


0.38 


0.96 


0.60 


2.43 


Y 


Y 


Y 


15 


M6 


1.74 


0.00 


0.01 


0.01 


0.01 


0.03 


N 


Y 


Y 


16 


K4: 


1.83 


0.32 


0.09 


0.44 


0.12 


0.61 


N 


Y 


Y 


18NE 


M5 


0.60 


0.26 


0.15 


0.45 


0.25 


0.83 


Y 


Y 


Y 


18SW 


M2: 


0.76 


0.26 


0.12 


0.41 


0.27 


0.83 


W 


Y 


Y 


19 


MS: 


1.36 


0.25 


0.0 


0.25 


0.00 


0.25 


Y 


Y? 


Y 


20 


M2 


0.48 


0.38 


0.04 


0.43 


0.13 


0.63 


Y 


Y 


Y 


2NE 


MS 


0.20 


0.25 


0.00 


0.25 


0.18 


0.48 


W? 


Y? 


Y? 


2SW 


M5 


0.20 


0.25 


0.12 


0.40 


0.19 


0.67 


Y 


Y? 


Y? 


1 


Ml 


1.30 


0.25 


0.00 


0.25 


0.00 


0.25 


N 


Y 


7 



^ Spectral types flagged with colons are derived from somewhat noisy spectra, resulting in less certain 
derived spectral types. 

Spectra are extremely noisy, and lack data in the J band; measurements of vj and A{J-H) arc 
impossible, and other measurements are very low quality. 



Tabic [5] and Figs. [6]and[7l 

3.4. Calculating T^ff and L^oi 

To compare the properties of B59's YSOs with theoret- 
ical models of pre-main sequence evoluti on, we adopted 
the Sp T to Te// conversion developed by iLuhman et al.l 
()1998[ ). The Te// implied by each source's spectral type 
is shown in Table |6l In the spectral type range occupied 
by most B59 members, our ±2 subclass spectral type 
uncertainty corresponds to an ^ 300 K uncertainty in 
temperature. 

With extinction and veiling estimates for each candi- 
date B59 member, we are able to calculate stellar lu- 
minosities f r om th eir NIR photometry. We followed 
IHillenbrandl (|1997f) in calculating each YSO's bolomet- 
ric luminosity as: 



Log{Ltot/LQ) = 0.4 x 
f4.75 -{H -Ah - DM + BCh 



(6) 



where DM is t he distance modulous to B59 (d=130 pc, 
DM = 5.569; iLombardi era!] [20061) and BCh is the H 
band bolometri c correction for each s t ar's s pectral type 
as tabulated bv lKenvon fc HartmannI (|1995[ ): the values 
adopted for each star are included for reference in Table 

m 

We also calculate the purely stellar component of each 
YSOs bolometric luminosity by removing the H band 
veiling flux, expressed in magnitudes (see ^3.21 for the 
transformation between the veiling flux ratio in a given 
band (e.g, th) and the photometric excess in that band 
(e.g., A H): 



(4.75 - {H + AH - Ah - DM + BCh, 



(7) 



Log{Lstei/LQ) = 0.4 x 



We adopt the measured 2MASS H band magnitude for 
isolated B59 members in this calculation, enabling us to 
perform an identical analysis of YSOs with known spec- 
tral types in other well studied star forming regions. Ob- 
jects with nearby (<4") companions may be blended in 
2MASS due to that survey's relatively large pixels (^2"); 
for those sources ([BHB2007] 2, 8 NW & SE, 13 NW, 
18 NE & SW; [CLR2010] 2 NE & SW), we adopt the 
synthetic H band magnitude measured from our spectra, 
noting that the SpeXTool reduction pipeline produces 
spectra with a first-order correction for slit losses derived 
from the observations of the telluric standard. 

Both luminosity estimates (total and stellar-only) are 
listed in Table |6] and shown in Figure [51 Accounting for 
the uncertainties in each of the input parameters results 
in a total uncertainty for these luminosity estimates of 
0.3 dex. The difference between the total and stellar 
luminosity estimated for a given YSO is usually less than 
0.2 dex, indicating that the veiling correction, while not 
negligble, does not exceed the overall error in each YSO's 
location in the HR diagram. 

3.5. Inferring B59's Median Stellar Age 

We have used the pre- main sequence (PMS) evo l ution - 
ary models calculated bv lD'Antona fc Mazzitelhl ()1994( ) 
and up dated in 1998 (DM98 hereafter), and those calcu- 
lated bv lBaraffe et al.l()1998[ ). to infer the age and mass of 
each candidate B59 YSO. These parameters were deter- 
mined by interpolating between model points with tem- 
peratures and luminosities similar to that of each YSO. 
These theoretical models do not account for accretion lu- 
minosity, so for this comparison we adopted the veiling 
subtracted value of Lstei for each YSO. Table [6| presents 
masses and ages derived for candidate B59 members with 
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TABLE 6 

Inferred Stellar Parameters 



Source 



'-iff 



BCh Log 



Lta 



Log 



DM98 Age 

(Myr) 



DM98 Mass 

(M©) 



B98 Age 

(Myr) 



B98 Mass 
(Mq) 



BHB2007 


1 


4060. 


2.11 


0.02 


-0.10 


0.79 


0.49 


4.37 


1.14 


BHB2007 


2 


3346. 


2.37 


-0.03 


-0.03 


<0.07 


0.20 


<1.0 


0.60 


BHB2007 


3 


2841. 


2.58 


-1.38 


-1.38 


2.59 


0.09 


<1.0 


0.07 


BHB2007 


4 


3009. 


2.51 


-0.75 


-0.89 


1.04 


0.14 


<1.0 


0.14 


BHB2007 


6 


3514. 


2.27 


-0.14 


-0.14 


0.15 


0.24 


<1.0 


0.62 


BHB2007 


7 


4350. 


2.02 


0.12 


-0.22 


2.71 


0.75 


12.02 


1.16 


BHB2007 


9 


4350. 


2.02 


-0.37 


-0.67 


23.4 


0.79 


50.11 


0.77 


BHB2007 


12 


3009. 


2.51 


-1.12 


-1.27 


2.77 


0.14 


1.38 


0.1 


BHB2007 


13NW 


3514. 


2.27 


0.03 


-0.03 


0.088 


0.23 


<1.0 


0.68 


BHB2007 


14 


4350. 


2.02 


0.05 


-0.24 


3.03 


0.77 


12.59 


1.13 


BHB2007 


15 


2841. 


2.58 


-0.90 


-0.91 


0.28 


0.11 


<1.0 


0.12 


BHB2007 


16 


4590. 


1.98 


0.07 


-0.09 


2.77 


0.88 


13.18 


1.24 


BHB2007 


18NE 


3009. 


2.51 


-0.78 


-0.94 


1.28 


0.14 


<1.0 


0.13 


BHB2007 


18SW 


3514. 


2.27 


-0.89 


-1.04 


6.71 


0.34 


13.2 


0.46 


BHB2007 


19 


3346. 


2.37 


-1.37 


-1.47 


16.6 


0.26 


19.1 


0.27 


BHB2007 


20 


3514. 


2.27 


-0.44 


-0.60 


1.37 


0.30 


3.47 


0.50 


V359 Oph 


4060. 


2.11 


-0.36 


-0.36 


2.42 


0.61 


10.47 


1.02 


Lk Ha 34E 




4060. 


2.11 


-0.19 


-0.39 


2.83 


0.63 


10.96 


1.00 



CLR2010 
CLR2010' 



2NE 
2SW 



3346. 
3009. 



2.37 
2.51 



-0.59 
-0.70 



-0.72 
-0.84 



1.28 
0.54 



0.22 
0.13 



2.19 
<1.0 



0.33 
0.14 



reliable spectroscopic and photometric information. 

The stellar ages presented in Table [6] allow us to iden- 
tify the characteristic timescale over which B59 has been 
actively forming stars. We calculate this timescale as the 
median age of all likely B59 members with reliable age es- 
timates (i.e., ah YSOs hsted in Table H except [CLR2010] 
2 NE & SW, which may not be linked to B59). We prefer 
the median to a standard mean because it is more robust 
to outliers with extreme ages (e.g., [BHB2007] 19). B59's 
median stellar age, as inferred from DM98 PMS tracks 
using veiling-corrected luminosity estimates {Lstei), is 2.6 
Myrs. 

Our confidence in the age estimates we calculate here is 
reinforced by our abilit y to independently repro duce the 
age estimate derived bv lLuhman fc Riekg (119991) for their 
sampl e of p Oph members. Following iLuhman fc Riekd 
(jl9 99f) by comparing each star's total (non-veiling cor- 
rected) luminosity to the DM98 tracks, and adopting 
their estimated distance of 166 pc. (to — M = 6.1), 
we measure a median stellar age for p Oph of 0.21 
Myrs, consistent with t he '^O.S Myr median stellar age 
ILuhman fc Riekd (|1999fl derived. 

We also measured median stellar ages for other, well 
studied young clusters to provide an indication of B59's 
relative, rather than absolute, age. These cluster age 
measurements require YSOs catalogs with measured 
spectral types and JHK photometry to estimate ex- 
tinction and veiling corrected stellar luminosities using 
the same algorithms we applied to B59's YSOs. Stud- 
ies by Luhnian and collaborators provide such catalogs 
for YSOs in Taurus, R h o Ophiuchus, and Chameleon 
(ILuhman fc Riekd Il999t ILuhman et alj l2006t iLuhmanI 
I2007D . We have placed these objects on the HR diagram 
using the same routines applied to our B59 data, adopt- 
ing the following distances to each r egion: Taurus - 140 
PC Oph - 126 PC Cha - 170 pc (jKenvon et all 120081 : 
iWilking et aIll2008l:lLu"hman 2008). Table [7] presents the 
me dian stellar age impli ed for each cluster by the DM98 
and lBaraffe et alJ () 19981 ) PMS tracks, and Figure[9]shows 
the locations of each region's YSOs in the HR diagram. 

Our analysis of the ages of YSOs in other regions. 



however, lacks the iterative refinement of each YSO's 
veiling estimate according to the agreement between the 
strengths of spectral features in the target and artificially 
veiled standard spectrum. Instead, for YSOs in other re- 
gions, we simply adopt the mean CTTS J band veiling 
measured by Cieza et al. (2005), which implicitly as- 
sumes that CTTSs dominate the population of each star 
forming region. This J band excess is an overestimate of 
the minimal veiling flux detected from WTTSs, however, 
such that we will under-estimate WTTSs intrinsic stel- 
lar luminosity, and thus overestimate their age. Our age 
estimates for Taurus, p Oph and Chameleon will there- 
fore be slightly biased towards older ages in proportion 
to the number of WTTSs in the sample: we expect this 
to be a minimal effect for Taurus and p Oph, but poten- 
tially of some importance for the age we measure for the 
Chameleon star forming region. 

3.6. Uncertainties in Derived Ages 

Our measurement of B59's median stellar age is subject 
to a few distinct statistical and systematic uncertainties. 
These affect the relative and absolute precision of our 
measurement in different ways, so we discuss them each 
below. 

3.6.1. Random/ Statistical Age Uncertainties 

Our measurement of cluster median stellar ages are 
susceptible to (presumably) random uncertainties due to 
observational error and sampling effects. Age uncertain- 
ties introduced by observational errors are straightfor- 
ward: errors in a YSO's temperature and luminosity will 
shift its location in the HR diagram, corresponding to 
an error in the inferred mass and/or age. In the tem- 
perature/luminosity regime occupied by B59's YSOs, lu- 
minosity is primarily sensitive to a star's age, while its 
temperature is a better indicator of its mass. We there- 
fore estimate the age errors that could be introduced by 
uncertainties in our luminosity estimates: we find that a 
0.3 dex luminosity error corresponds to an age error of 
~125%. This is, however, the error in any individual age 
estimate: as we are determining the median stellar age of 
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TABLE 7 
Median Cluster Ages 



Cluster 


Nstars 


DM98 Median 


B98 Median 


Statistical 


Systematic 






Age (Myrs) 


Age (Myrs) 


Uncertainty (Icr) 


Uncertainty 


B59 


18 


2.59 


4.37 


29% 


100-175% 


Taurus 


224 


1.53 


2.63 


9% 


100-175% 


p Oph 


61 


1.04 


3.63 


16% 


100-175% 


Chameleon 


185 


2.59 


3.02 


9% 


100-175% 



o 

+ 




("h'jV HD 181596 
' ' ' K5 11I,A,=6.5 



\i)^ HD 196610 

M6 III. A^=3.25 

iTVy ICLR2010]3 



^ P^rvS :::::: 

iHW MYCep 







Pfi.j. G1686 






1.0 



2.5 



1.5 2.0 

Wavelength (microns) 

Fig. 5.— Spectroscopically identified giant stars in B59 
and the Pipe Nebula. Top: K & M giants and super- 
giants - [BHB2007] 5, [CLR2010] 3, 7, & 9, shown in 
comparison to spectra of HD 181596 (K5 III), HD 156014 
(M6 HI), a nd MY Cep (M7 I) from the IRTF spec- 
tral library (jCushing et al.l 120051 ) . HD 181596 and HD 
196610 have been artificially reddened to better match 
the target spectra. Middle: Carbon stars - [CLR2010] 
8 show n for comparision with the reddened spectrum of 
R Lep (jLanpon fc Woodl[2000l ). Bottom: OH/IR stars 
- [CLR2010] 4, 5 & 6, as well as [BHB2007] 17, shown 
for c omparison with the spec trum of OH/IR star Gl 686, 
from lLanconfc WoodI (poM ). 

all the YSOs in ea ch cluster , we can expect that this er- 
ror will scale like \/N duster- We calculate an estimate of 
the statistical uncertainty in each cluster's median stellar 
age as /I? , and present the resulting uncertainty in 

Table m 



Sampling errors are perhaps the more subtle effect, but 
could be important for a sparse cluster like B59. A sim- 
ple demonstration of this effect can be seen by examining 
the ages inferred for the 18 B59 YSOs listed in Table H 
The median age of these B59 YSOs, as inferred from the 
DM98 models, is bounded by the 2.42 Myr and 2.59 Myr 
ages estimated for V359 Oph and [BHB2007] 3, respec- 
tively the 9*'* and 10*'' oldest YSOs in the fist. The rela- 
tively narrow age range bounding the median, however, 
is somewhat misleading. There are at least three B59 
YSOs that were too heavily embedded to be observation- 
ally accessible to SpeX: [BHB2007] 10 and 11, classified 
as Class O/I and Class I sources bv lBrooke et a l. (2007), 
and 2M17 112255, discovered and classified as a Class I 
system bv iRiaz et al.l ()2009[ ). These YSOs are presum- 
ably quite young: arbitrarily assigning them ages of 0.5 
Myrs would shift the age ordering such that the YSO 
defining B59's median stellar age would be [BHB2007] 
20, with an age of 1.37 Myrs according to the DM98 
models. This would represent a ^50% age error, some- 
what lar ger than the 30% age estimate predicted by the 
simple \/Nciuster calculation described above. We lack, 
however, a robust statistical description of the effect, so 
merely caution the reader that the formal statistical un- 
certainties reported in Table [7] may be underestimates 
for a sparse cluster like B59. 

3.6.2. Systematic Age Uncertainties: Veiling Corrected 
Luminosities 

One systematic effect which could affect our age es- 
timates concerns the treatment of veiling flux in calcu- 
lating each YSO's luminosity. Deriving ages from each 
YSO's total luminosity, rather than the veiling corrected 
stellar luminosity, typically reduces the YSO's implied 
age by ^1.5 Myrs over the range of ages considered here. 
The median stellar age implied for B59 by comparing 
each YSO's total, non- veiling corrected luminosity esti- 
mate to the DM98 tracks is 0.87 Myrs, as opposed to 
the veiling corrected age estimate of 2.6 Myrs. This 
represents a nearly 75% decrease in the cluster's me- 
dian stellar age, demonstrating that correctly accounting 
for veiling luminosity is necessary to accurately infer a 
YSO's age. As the PMS tracks considered here do not 
account for emission from accretion or a circumstellar 
disks, we believe the correct approach is to infer ages 
from veiling-corrected luminosity estimates. We caution 
the reader, however, that many other investigators adopt 
non-veiling corrected bolometric luminosity estimates to 
derive YSO ages: as demonstrated above, this difference 
in approach results in age differences of 50-100%, and 
should be carefully considered in reconciling our results 
with other studies in the literature. 

3.6.3. Systematic Age Uncertainties: Adopted PMS Tracks 
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Fig. 6. — H and K band spectra ol confirmed or candidate B59 nienibers.^BHB2007] 8NW & SE are not shown, as 
their spectra are very low signal-to-noise. 
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Fig. 7. — J band spectra of confirmed or candidate B59 
members with visible J band emission lines. 



As shown by iHillenbrand et ahl (|2008[ ) , the choice of 
PMS models from which to infer YSO ages has a signif- 
icant affect on the absolute age derived for a given clus- 
ter. [Hillenbrand et al. (200 8) found that the PM S mod- 
els computed by DM98 and lBaraffe et al.l (| 19981 ) predict 
the youngest and oldest ages, respectively, for a YSO of 
a given temperature and luminosity. We therefore use 
these two sets of model tracks to span the full range of 
B59's potential absolute age. Using the same veiling cor- 
rected luminosities derived above for B59's YSOs, but 
inferring ages from the a = 1 PMS tracks calculated by 
iBaraffe et al.l ()1998[) . results in a median stellar age for 
B59 of 4.4 Myrs, or roughly 70% older than the median 
stellar age inferred from the DM98 models. 

Model dependent age differences are often attributed 
to luminosity offsets between equal age isochrones in two 
sets of PMS tracks. This suggests that, while the abso- 
lute ages predicted by two sets of PMS tracks may differ, 
they should produce similar relative age orderings for a 
given sample of clusters. Table [7] reveals, however, that 
this simple expectation is not necessarily borne out in 



practice: both sets of models indicate that B59 has the 
oldest median stellar age, but the DM98 models imply 
tha t p Oph is younger th an Taurus and Chameleon, while 
the IBaraffe et al.l (|1998[ ) models identify p Oph as older 
than Chameleon and Taurus. 

These different relative age orderings arise from differ- 
ences in the stellar population sampled in each cluster, 
combined with distinct isochrone slopes predicted by the 
two sets of PMS tracks. The temperature dependent off- 
set between the 5 Myr DM98 and Baraffe isochrones is 
shown in the left panel of Figure [TOl The Baraffe models 
imply consistently older ages for a given HR diagram lo- 
cation, but the age difference increases rapidly towards 
higher Tg//. This means that the relative age orderings 
of single stellar populations may not be preserved across 
models if the YSOs observed sample different tempera- 
ture ranges. This is precisely the cause of the different 
relative ages inferred for p Oph in the DM98 and Baraffe 
models: as demonstrated in the right panel of Figure [TUl 
the median temperature of the p Oph cluster stars is sig- 
nificantly larger than that of the other clusters analyzed 
here. 

3.6.4. Systematic Age Uncertainties: Cluster Distances 

Uncertainties in a cluster's distance can also introduce 
errors in the inferred age: a closer distance implies a 
lower luminosity, and thus an older age, for a YSO of 
a given apparent magnitude. A concrete example of of 
this effect is provided by p Oph, whose distance has re- 
cently been revised to ~125 parsecs from a previous esti- 
mate of ~^160 parsecs. The 0.3 Myr es timate of p Oph's 
age derived bv lLuhman fc Riekd (|1999[ ) was based on the 
160 parsec distance estimate. Replicating their analysis 
with a revised distance of 125 parsecs increases the de- 
rived median stellar age to 0.63 Myrs, twice as old as 
the previous estimate. Estimates of the distance to B59 
and/or the Pipe Nebula are scarce: the two most recent 
distances measured for the Pipe are 130^20 parsecs and 
145ibl6 parsecs, obtain ed bv lLombardi et all (|2006D and 
lAlves fc Francol (I2007f). re spectively We have adopted 
the ILombardi et al.l (|2006t ) distance for the luminosity 
an d age measurements pre sented here, but we note that 
the lAlves fc Francol ()2007f ) distance implies in a median 
stellar age (inferred from veiling corrected luminosities 
and DM98 tracks) for B59 of 1.7 Myrs, or decreased by 
roughly 35% from our primary age estimate. 

3.6.5. Total Error Budget for B59's Age Estimate 

In summary, our measurement of B59's absolute age 
is subject to statistical errors at the ^30% level, such 
that B59's median stellar age (adopting veiling corrected 
luminosities and inferred from the DM98 models) is 
2.6±0.8 Myrs. Several systematic uncertainties could 
shift the median stellar age derived for B59, however, by 
as much as 100%. Adding these effects in quadrature, we 
find that systematic uncertainties could introduce errors 
as large as 150% into our measurement of B59's absolute 
age. Including potential systematic effects therefore re- 
vises the error budget for B59's median (DM98) stellar 

age to 2.6^2:6 ^yrs. 

4. DISCUSSION 

4.1. The Star Formation Timescale in B59 
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Fig. 8. — B59 HR diagram calculated using JHK photometry. Likely B59 members (including previously studied 
members V359 Oph & Lk Ha 345) are shown twice; red dots indicate their total bolometric luminosity, and black dots 
indicate their stellar luminosity after subtracting H band veiling. Objects classified as background giants are shown 
with open diamonds, as are [BHB2009] 8SE & NW, whose spectra are too noisy for reliable HR diagram placement. 
The DM98 model grid is shown for comparison as solid and dotted lines. 



The prompt star formation paradigm agrees well with 
the evolutionary timescales expected for dense cores that 
lack significant internal or external support. Acted on by 
gravity alone, a core should collapse in roughly one free- 
fall time, Tff = (37r/32G/5)^/^. Hydrodynamic models of 
undriven turbulent molecular clouds predict only slightly 
longer evolutionary timescales, as the t urbulent sup- 
port di ssipates on ~l-2 free -fall timescales (jKlessen et al.l 
I2000f ). lOffner et all (|2008[ ). for example, find that star 
formation activity ceases after 0.75 t// in their simula- 
tions of star formation within undriven turbulent molec- 
ular cores. 

The uncertainties associated with our measurement of 
B59's age are sizable, but the most straightforward in- 
terpretation of this result is that star formation activity 
began in B59 nearly 3 Myrs ago. The presence of heav- 
ily embedde d (and presu mably young) Class O/ I objects 
within B59 ([Brooke et a l. 2007; Riaz et al.ll2009D also im- 
plies that this star formation activity has persisted to the 
present day. If the oldest stars in B59 formed from the 
dense core currently associated with the cluster, these 
observations would indicate a lower limit of 3 Myrs for 
the core's lifetime. A star formation timescale for B59 
comparable to 3 Myrs is by no means startling in ab- 



solute terms. Indeed, while the stellar content of B59 
appears somewhat older than that of the p Oph, Taurus, 
and Chameleon star forming regions, the inferred star 
formation timescales agree to within a factor of two. 

Expressing B59's age in units of the core free-fall time 
and the core's crossing time is also useful for interpret- 
ing the core's physical evolution. Adopting the den- 
sity (n(H2) = 10'' cm-3; p = 3.5 x 10"^" g cm'^), 
radius (R = 0.19 parsecs — 5.9x10^^ km), and veloc- 
ity dispersion (cti,=0.37 km/sec, calculated from 0.86 
km/sec FWHM) tabul ated ifor the B59 dense core by 
iRathborne et al.l ()2009D . we calculate Tff ~0.35 Myrs 
and Tcross ~0.5 Myrs. These calculations describe B59's 
current state; if B59 is evolving from a lower density to 
higher density configuration, these timescales would have 
been somewhat longer in the past than they are today. 
Intensive observations of the dense gas within B59 have 
not revealed any kinem atic signatures of rapid collapse 
(jRathborne et a l. 20091. If B59 is evolving relatively qui- 
escently, the timescales we derive today are likely com- 
parable to those of B59's recent past. 

The uncertainties on B59's age permit that star forma- 
tion may have begun as recently as the core's last dynam- 
ical time, or as long as ^13 dynamical timescales ago. 
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The maximum likelihood value, however, suggests that 
B59 has been forming stars for f«6 crossing or free fall 
times. Extinction maps of B59 also indicate that the core 
is spatially coherent, having not undergone significant 
fragmentation as it formed stars ([Roman- Zuhiga et alj 
[2009) . Taken at face value, the maximum likelihood 
value of B59's star formation timescale suggests that 
B59 could join a short, but growing, list of dense cores 
that appear to have resisted global collapse and sus- 
tained ongoing star formation for several dynamical 
times. [Swift & Welch (2008) examined the dense gas 
and stellar population of the L1551 star forming region, 
demonstrating that L1551 has been forming stars for ~5 
Myrs, or >5 r//; they identify stellar feedback as a likely 
source of support against collapse within this region. As 
we find for B59, |Swift fc Wel cb[ ((2008I ) also find that 
L1551 contains stars over a broad range of evolutionary 
stages (Class through III). Similar re sults have been 
report ed for the young cluster AGFL961: lWilliams et al.l 
([2009( ) conclude that either disk evolution proceeds more 
rapidly in clustered environments, or that star formation 
persists in cluster s over several dynamical timescales. 
Finally, I Tan et ahl ([20060 present observationally-based 
calculations that suggest ages of ~3-4 r// for several rich 
young clusters that are still forming stars (ie, p Oph, IC 
348, ONC, etc.). These measurements suggest that some 
dense cores can sustain star formation over a few dynam- 
ical timescales, including in relatively low-mass cores like 
B59, which does not appear to host massive stars. 

What physical mechanism could be supporting these 
cores against collapse, extending their lifetimes beyond 
that expected for cores in free fall coll apse? As noted 
above, one possibility is stellar outflows: iMatzneii ([2007f ) 
calculate that outflows may support clusters against col- 
lapse for several crossing times. This mechanism could 
explain B59' s relat ive s tability, pa r ticula rly given that 
lOnishi et all ([19990 and iRiaz et al.l ([2OO90 find evidence 
for an active outflow within the B59 core. A recent 
ammonia survey of the Pipe, however, finds little evi- 
dence for significant non-thermal motions within B59's 
dense gas, as would be expected if bulk outflow mo- 
tions are the core's main supp o rt ag ai nst gravitational 
collap se ([Roman-Zuniga et al.l 120091 : iRathborne et al.l 
120080 . High-resolution observations of the kinematics of 
B59's dense gas are required to clarify the extent to which 
it (and, by extension, outflows from previous generations 
of protostars) could support B59 against collapse. 

The collapse of a molecular core may also be strongly 
effected by radiative feedback from the young protostars 
it contains. Krumholz, Klein & McKee (2007) used sim- 
ulations of massive molecular cores to demonstrate that 
protostars can contribute significant heating to a collaps- 
ing core, thereby suppressing fragmentation. This heat- 
ing does not, however, appear to slow mass accretion and 
collapse within the core: much of the radiation 'leaks' 
out of optically thin holes in the core (Krumholz et al. 
2009). Offner et al. (2009) extended this investigation 
into the regime of low- mass star formation, finding a sim- 
ilar reduction in stellar multiplicity due to reduced frag- 
mentation, as well as a significant reduction in the mass 
accretion rate of the earliest formed protostar. These 
simulations suggest that radiative feedback could play 
some role in extending B59's star formation timescale. 

Magnetic fields are a final physical mechanism that 



is often invoked to support cloud cores against col- 
lapse. Magnetic fields have been a fundament al com- 
ponent of star formation models for dec ades ("N akanol 
[l98l IShu etHI [19871 iMouschovias fc Cio lck 1999), but 
observational constraints on the strengths of magnetic 
fields within molecul ar clouds are difficult to obtain (e.g., 
iCrutcher et al.|[20090 . Past measurements suggest mag- 
netic fields may contrib ute appreciable amounts o f sup- 
port to molecular cores ([Troland fc Crutcheill20080 . and 
recent polarization measurements in the Pipe Nebula in- 
dicate the presence of a ^1 7 p,G magnetic field within 
B59 itself ([Alves et al.ll2OO80 . This suggests that mag- 
netic fields may be providing B59 with some support 
against collapse, but this is a relatively weak statement: 
magnetic fields can extend a dense core's lifetime arbi- 
trarily, depending on the field strength adopted. 

4.2. The Star Formation Efficiency of B59 

As the discussion above demonstrates, there are several 
mechanisms which could plausibly explain B59's maxi- 
mum likelihood star formation timescale. To provide ad- 
ditional leverage for comparing observations of B59 with 
predictions of theoretical models of star forming cores, 
we now consider B59's star formation efficiency in addi- 
tion to its characteristic timescale. 

B59's present-day star formation efficiency (SFE) can 
be characterized as the mass of the cluster's stellar pop- 
ulation divided by the core's initial gas mass. The to- 
tal mass of all B59 YSOs analyzed here is '-^S M© ac- 
cording to the DM98 models, and 12 M(t, accor d ing to 
the Baraffc et al. ( 1998) models. iBrooke et all ([2OQ70 
and iRiaz et al.i (|2009il also provide preliminary mass es- 
timates for [BHB2007] 10, 11, and 2M12112255, the B59 
YSOs that were too heavily embedded to be observed 
with SpeX. These embedded YSOs contribute an addi- 
tional 0.7 M© to B59's total stellar mass, bringing B59's 
total stellar mass to ~9-13 Mq. 

B59's initial gas mass can be approximated as the sum 
of its stellar and gas masses. Extinction maps indicate 
that B59 contains ^20 Mq of dense gas, implying an 
initial core mass of ^29-33 Mq. This suggests B59's 
present-day SFE is ^30-40%, depending on the choice 
of models used to infer stellar masses. Combining this 
SFE with B59's median stellar age produces an estimate 
of the SFE/t// of the B59 dense core. B59's sizable 
age uncertainty translates into a large range of allowed 
SFE/t// values: 3-40%, though the probability distri- 
bution is strongly asymmetric, peaking at a maximum 
likelihood value of ~6%. 

The SFE/Tff values predicted by theoretical mod- 
els for star forming cores cover a large range, with 
slower star formation rates corresponding to models in- 
corporating some additional physical mechanism (driven 
turbulence, magnetic fields, outfiows, etc.) to sup- 
port cores against collapse. Unsupported models typi- 
cally produce SFE/r// values ~ 10-20%, while models 
that incorporate additional support often predict pro- 
portionally lower SFE/ryj values. As one example, 
the moderately supercritical (i e, prone to gravitational 
collap se) magnetic models of I Vazquez- Semadeni et al.l 
([2OO50 produced SFE/r// values ~ 5%, whereas ear- 
lier non-magnetic model s by the same authors produce d 
SFE/t// values ~15% ([Vazauez -Scmadcni et al."200aV 
The magnetically subcritical models of ; Nakamura fc Lil 
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(|2008D , which also incorporate feedback from stellar out- 
flows, produce SFE/r// values of 1% or less^ 

The theoretical models calculated by iPrice &: Bate! 
(|2009t ) are arguably the best match to B59's actual phys- 
ical structure. These models (indicated hereafter as 'the 
PB09 model grid') incorporate both magnetic fields and 
radiative feedback, and follow the evolution of a cloud 
core whose mass (50 Mq) and radius (~0.2 pc) are rea- 
sonably close to the values which likely characterized B59 
prior to the onset of active star formation. The PB09 
model grid spans magnetic fields of various strengths, 
but all models were supercritical, with mass-to-flux ra- 
tios ranging from 3 to c». Though the models were only 
allowed to evolve over a timespan corresponding to 1.5 
Tff, or 0.35-0.5 r// following the formation of the first 
star in the core, the implied SFE/t/^ values range from 
^-^30% for the model lacking both magnetic fields and 
radiative feedback, to ~ 6% for the model incorporat- 
ing radiative feedback and the strongest magnetic field 
(Bo=65/iG). In analyzing the range of star formation ef- 
iiciencies produced by their model grid, PB09 note that 
magnetic fields appear to play a more important role 
than radiative feedback in quenching star formation. 

B59's maximum likelihood SFE/r^j value of ^-^6% is 
contained within the range of values predicted by the 
PB09 models, but B59's observed magnetic field strength 
(B^17 /iG) is only ^1/4 that adopted in the least effi- 
cient PB09 model (B-65 fiG). The PB09 models and the 
present-day B59, however, represent significantly differ- 
ent time periods in the evolution of a star forming core. 
It would be interesting to see what properties the PB09 
models might have if they were allowed to evolve to ^5 
Tff. would the magnetic field diffuse out of the core, such 
that the measured magnetic field strength would better 
agree with what we observe in B59 today? 

It is too early to draw detailed conclusions, but the 
agreement between B59's observed properties and the 
predictions of the PB09 model grid are intriguing, and 
suggest that magnetic support and radiative feedback 
help support B59 against collapse. 

4.3. The Spatial Extent of Star Formation in Pipe 

The suggestion that magnetic fields may be important 
for understanding star formation within the Pipe Nebula 
is also consistent with the results of t he polarizati o n sur - 
vey conducted bv lAlves et all (|2008l ). lAlves et al.l (|2008f ) 
measure a significant spatial gradient in the component 
of the Pipe's magnetic field aligned with the plane of the 
sky, rangin g from ^17 /xG in B59 to ~65/xG in the Pipe's 
'Bowl'. As lAlves et al.l()2008[ ) indicate, this gradient cor- 
relates with star formation activity: the weakest, most 
poorly aligned field vectors lie in B59, where the bulk of 
the Pipe's star formation activity is taking place, while 
areas with stronger, better aligned magnetic field vectors 
appear devoid of star formation activity. 

In this context, it is interesting to note the loca- 
tion of [CLR2 10] 2 the candidate YSO identified by 
iForbrich et~all (|2009 ") in their MIPS survey of the Pipe, 
and subsequently confirmed as a likely YSO by our spec- 
troscopic analysis. This candidate YSO lies in the Pipe's 
'Stem', where the magnetic field appears to be stronger 
than in B59, but not as strong as in the Bowl. This YSO 
also displays a rich emission line spectrum, suggesting it 
is a young, heavily accreting YSO, which would be con- 



sistent with a picture in which star formation activity was 
just beginning in the Neck of the Pipe. We note, however, 
that [CLR2010] 2 does not appear to be heavily embed- 
ded; indeed, [BIIB2007] 10 and 11 appear to be in sig- 
nificantly earlier evolutionary stages than [CLR2010] 2, 
making it difficult to conclude that the onset of star for- 
mation activity is simply progressing linearly from B59 
towards the Bowl of the Pipe. 

5. SUMMARY 

1. We obtained and analyzed JHK SpeX spectra for 
20 candidate B59 YSOs, and 10 additional mid-IR 
bright sources throughout the larger Pipe Nebula. 
We measured temperature and luminosity sensitive 
spectral indices from these spectra, and determined 
the extinction, NIR color excesses, and spectral 
type required to self-consistently reproduce each 
target's spectral energy distribution. This anal- 
ysis leads us to re-classify [BHB2007] 5 and 17, 
previously identified as candidate B59 members, 
as background giants, and to identify [CLR2010] 
2 SW as a promising candidate YSO in the stem 
of the Pipe. None of the mid-IR bright candidates 
outside of the Stem of the Pipe were found to be 
YSOs 

2. For each candidate B59 YSO with reliable spectra 
and photometry, we derived effective temperatures 
and luminosities, with and without correction for 
non-photospheric flux. Using these measurements 
to place each YSO in the HR diagram, we inferred 
ages and masses for each YSO from theoretical 
pre-main sequence evolutionary tracks. B59's me- 
dian stellar age, as inferred from the DM98 mod- 
els and accounting only for statistical uncertainties, 
is 2.6±0.8 Myrs. Including potential systematic 
effects revises the error budget for B59's median 
(DM98) stellar age to 2.6^2:6 ^yrs. 

3. Comparing the median stellar ages implied for sev- 
eral star forming regions by different pre-main se- 
quence evolutionary tracks demonstrates that the 
construction of even a relative age ordering requires 
careful attention to sampling of each region's stel- 
lar population. Despite this sensitivity, both sets of 
tracks studied here imply that B59 is at least as old 
as, if not slightly older than, the more well-studied 
p Oph, Taurus, and Ghameleon star forming re- 
gions. The dominant source of error in this con- 
clusion is the uncertainty in the distance adopted 
to each cluster: we minimize other potential sys- 
tematic errors by applying a consistent set of algo- 
rithms to a homogeneous set of catalogs of young 
stars in each cluster. 

4. The maximum likelihood median stellar age we 
measure for B59, and the observed properties of 
the dense gas in the region, suggests that the B59 
dense core may have been stable against collapse 
and actively forming stars for roughly 6 dynami- 
cal timescales, with a star formation efficiency per 
dynamical time of ^-^6%, though the uncertainty in 
our age measurement permits a significant range of 
possible SFE/t// values (3-40%). These proper- 
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ties agree well with recent star formation simula- 
tions that incorporate sub-critical magnetic fields 
and radiative feedback from protostellar heating. 
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APPENDIX 
TWA SPECTRAL STANDARDS 



To supplement the spectral atlas presented by iCushing et al.l (|2005[ ) with additional pre-main sequence spectral 
standards, we obtained SpeX SXD spectra of several members of the TW Hya moving group. Observed during the 
same nights as our B59 candidates, these spectra share the same instrumental set-up and reduction process as our 
target spectra. The H and K band portion of these spectra are shown in Figure 111! ordered such that earlier type 
objects are found higher in the figure. Spectral features sampled by the indices presented in H3.2l are also labeled. 

As most B59 members are heavily reddened, the z and J band portion of their observed spectrum was typically quite 
noisy, and thus not useful for measuring spectral types. This region contains most of the HI Paschen series, the Ca II 
triplet and the He I 10830 A line, however, which are often seen in emission in young, actively accreting stars. The 
0.8-1.3 lira SpcX spectrum of TW Hya, an actively accreting CTTS, is shown in Figure [T^for comparison with Figure 
[71 which shows 0.8-1.3 /im spectra of candidate or confirmed B59 members with J band emission features. 

To provide the community with moderate resolution NIR pre-main sequen ce spectral standards , we have made these 
spectra publicly available at |http : //www . astro . Cornell . edu/-kcovey/CoveyB59SpeX . tar . gz[ 
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